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ABSTRACT 

We model the acquisition of spin by dark-matter halos in semi-analytic merger trees. We 
explore two different algorithms; one in which halo spin is acquired from the orbital angular mo- 
mentum of merging satellites, and another in which halo spin is gained via tidal torquing on shells 
of material while still in the linear regime. We find that both scenarios produce the characteristic 
spin distribution of halos found in N-body simulations, namely, a log-normal distribution with 
mean ss 0.04 and standard deviation « 0.5 in the log. A perfect match requires fine-tuning of two 
free parameters. Both algorithms also reproduce the general insensitivity of the spin distribution 
to halo mass, redshift and cosmology seen in N-body simulations. The spin distribution can be 
made strictly constant by physically motivated scalings of the free parameters. In addition, both 
schemes predict that halos which have had recent major mergers have systematically larger spin 
values. These algorithms can be implemented within semi-analytic models of galaxy formation 
based on merger trees. They yield detailed predictions of galaxy properties that strongly depend 
on angular momentum (such as size and surface brightness) as a function of merger history and 
environment. 



Subject headings: galaxies:formation-galaxies:spiral 

1. Introduction 

Understanding the origin of spin in rotationally 
supported disk galaxies is clearly a crucial part 
of any theory of galaxy formation (Hoyle 1951). 
The general theory (Fall & Efstathiou 1980) re- 
produces galactic disks with roughly correct sizes, 
based on the assumptions that the gas originally 
had the same specific angular momentum as dark- 
matter (DM) halos do today, and that the infalling 
gas conserves specific angular momentum (Mestel 
1963). 

The origin of angular momentum in dark- 
matter halos can be understood in terms of lin- 
ear tidal torque theory in which protohalos are 
torqued by the surrounding shear field (Hoyle 
1951; Peebles 1969; Doroshkevich 1970; White 
1984; Barnes & Efstathiou 1987; Steinmetz & 



Bartelmann 1995; Porciani, Dekel, & Hoffman 
2001b). Many models for the formation of galac- 
tic disks have been proposed based on the general 
picture of Fall & Efstathiou. Some of these models 
construct disk galaxies solely from the properties 
of the final halos (Blumenthal et al. 1986; Flo- 
res et al. 1993; Dalcanton, Spergel, & Summers 
1997; Mo, Mao, & White 1998; van den Bosch 
2000), and others incorporate the mass-accretion 
histories of the halos (White & Frenk 1991; Kauff- 
mann, White, & Guiderdoni 1993; Cole et al. 1994, 
2000; Kauffmann 1996; Avila-Reese, Firmani, & 
Hernandez 1998; Firmani & Avila-Reese 2000; 
Somerville & Primack 1999; Somerville, Primack, 
& Faber 2001). 

The halo spin is a crucial element for modeling 
important physical quantities such as galactic-disk 
size, surface brightness, star formation rate, and 
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rotation velocity in semi- analytic models. How- 
ever, in none of the previous models is the angular 
momentum of DM halos assigned self-consistently 
based on their mass accretion and merger his- 
tories, primarily because an appropriate simple 
recipe for this buildup of spin has not yet been 
developed. The aim of this paper is to provide 
such a recipe. 

We explore two scenarios. In the first, halo 
spin is generated by the transfer of orbital an- 
gular momentum from satellites that merge with 
the halo (orbital-merger scenario). This scenario 
was suggested by Vitvitska et al. (2001), who also 
compare it with simulations in more detail than 
is done here. In the second, linear tidal-torque 
theory is applied to spherical shells of infalling 
matter (tidal-torque scenario). Each recipe is ap- 
plied in turn, within Monte Carlo realizations of 
halo merger histories generated using the extended 
Press-Schechter approximation and the method of 
Somerville & Kolatt (1999). We test and calibrate 
these scenarios by comparing their predictions for 
the distribution of spin with the known results of 
N-body simulations. 

The outline of the paper is as follows. Section 2 
presents our notation and defines the quantities of 
interest. Section 3 explores the orbital-merger sce- 
nario, while Section 4 addresses the tidal-torque 
picture. In § 5 we discuss the dependence of the 
spin distribution on mass, cosmology and redshift. 
In § 6 we study the dependence of the spin distri- 
bution on halo merger history. In § 7 we conclude 
and provide a brief discussion concerning the im- 
plications for disk sizes. 

2. Background 

The angular momentum of a halo, J, is com- 
monly expressed in terms of the dimensionless spin 
parameter (Peebles 1969) 

A = J^\E~\/GM 5/2 , (1) 

where E is the internal energy of the halo. In 
practice, the computation or measurement of this 
quantity, especially the energy, may be ambiguous. 
Furthermore, it introduces an undesired depen- 
dence on the specifics of the halo density profile, 
which depends on cosmology, redshift and halo 
mass (see Bullock et al. 2001) and may also de- 
pend on merging history (Wechsler et al. 2001b). 



Instead, following Bullock et al. (2000a), we will 
focus on the distribution of the variable 



V2MV vir R vil ' 

which is more straightforward to compute and 
does not explicitly depend on the density profile 
of the halo. For a singular isothermal sphere trun- 
cated at a radius R V1T , this parameter is equal to 
the traditional spin parameter, A' = A. They are 
also approximately equal for halos with a more re- 
alistic density profile: an NFW profile (Navarro 
et al. 1996) with concentration c v i r ~ 10. 

Many studies based on N-body simulations 
have revealed that the distribution of spin pa- 
rameter A is well fit by a log-normal distribution, 
and varies very little as a function of halo mass, 
redshift, or cosmological model. Bullock et al. 
(2000a) have demonstrated that the distribution 
of A' is also well fit by a log-normal distribution, 
namely, 

with X' Q = 0.035 and er A = 0.5 1 

3. The Orbital Merger Scenario 

In this section we explore the assumption that 
all of the orbital angular momentum of a merger 
is converted into the spin of the merger product. 
For simplicity we will consider all mergers as in- 
volving two bodies and refer to the more massive 
of the two progenitors as the halo and the other 
as the satellite. Following the merger tree algo- 
rithm (Somerville & Kolatt 1999), we neglect any 
mass or spin loss. We can safely neglect the in- 
ternal spin of the incoming satellite because even 
for equal mass mergers the spin is typically only 
<~ 10% of the orbital angular momentum. 

3.1. Encounter Parameters 

Consider a satellite whose center is at a distance 
r from the center of the halo and moving with a 

1 Note that in the log- normal distribution A is the spin cor- 
responding to the mean of ln(A'), and a\ is the standard 
deviation of ln(A'). We sometimes refer loosely to \' and 
c\ as the mean and standard deviation of the log-normal 
distribution. 
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velocity v relative to the halo. Then in center of 
mass coordinates, the energy and angular momen- 
tum are given by 
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1 2 M h m s 

Eorb = - G (4) 

2 r 

L rb = x r (5) 

where h and s denote the halo and satellite 
respectively and [i is the reduced mass, ^ = 
M h rn s / (Mh + m s ). We will assume that the satel- 
lite comes from outside of the halo and therefore r 
is at least as large as the virial radius of the halo 
(r > i? v ir)- For equal mass mergers (Mh = m s ) 
this implies that at some moment we can express 
the orbital energy in terms of the virial velocity of 
the halo 

Eorb = m *(r- Kir) , (6) 

which makes it clear that for a given energy there 
is a maximum orbital angular momentum L or b < 
fivRyit for the halo-satellite system. While this 
particular expression is only valid for equal mass 
mergers, the situation is similar for other encoun- 
ters. 

Binney & Tremaine (1987) introduce dimcn- 
sionless energy and angular momentum variables 
(Equation 7-85): 

E = J±±- ■ 1 = J^±_ ( 7 ) 

where r e and w rms are the half mass radius and 
rms velocity of the halo. For a halo with the 
NFW density profile and a concentration of « 15, 
r e = 0.3i? V i r and u rms = 1.2V^ir. For equal mass 
non-rotating spherical halos, the possible orbital 
parameters are shown in Figure 1. The dashed 
line excludes encounters that take longer then a 
Hubble time to merge. Adding the additional con- 
straint that the satellite starts entirely outside of 
the halo, (i.e. r > Rvh) leaves the hatched region 
as the permissible range of encounter parameters. 
Figure 1 is only valid for equal mass mergers; in 
general the possible values of E and L will depend 
on Mh and m s . 

Since we do not know the probability distribu- 
tions of E or L, we introduce a fudge factor / 
with which to parameterize our ignorance. Thus, 
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Fig. 1. — Initial encounter parameters. The pa- 
rameter space of orbit energy (E) and angular mo- 
mentum (L) for equal mass halos, following Bin- 
ney & Tremaine (1987). The region above the 
thick solid line is excluded by definition. The re- 
gion below the dashed line refers to encounters 
that will result in mergers within a Hubble time. 
The region below the thin solid line refers to ha- 
los that are separated by at least i? v ir- The cross 
marks the average value of L, and the diamond 
marks the value we use in our model. 

we will take the orbital angular momentum of the 
satellite-halo system to be given by 

Lorb = ffjy vil r. (8) 

If all permissible values for E or b are equally likely, 
we would deduce from Fig. 1 that a typical value 
for L is about 1, corresponding to / w 0.34. In 
principle, / could depend on mass and/or redshift 
but we will take it to be a constant and adjust it 
to get a good fit to the observed distribution of 
spin from N-body simulations; below we will find 
that / = 0.38 yields the best results. It would also 
be possible to determine the value of / empirically 
by analyzing N-body simulations (Vitvitska ct al. 
2001). 

We generate mass accretion histories based on 
the formalism of Somerville & Kolatt (1999) with 
the slight modification introduced by Bullock, 
Kravtsov, & Weinberg (2000b). We produce 500 
random realizations of the merger history of a 
halo, adopting a fixed value of /. The direction 
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Fig. 2. — The distribution of A' predicted by the 
orbital merger scenario for different values of / and 
7 m , compared to the simulation data of Bullock 
et al. (2000a) shown in the bottom-right panel. A 
choice of 7 m = 2. and / = 0.38 results in a A' dis- 
tribution identical to that seen in the simulations. 

of the orbital angular momentum, 9 (specifically, 
the angle between the orbital angular momentum 
vector of the satellite and the spin vector of the 
halo), is drawn at random from a uniform or Gaus- 
sian distribution. For the uniform case, we choose 
7 = cos((9) from a uniform distribution between -1 
and 1. For the Gaussian distribution, we choose 
values of 7 between -1 and 1 from: 

P( 7 )d 7 ex exp (- (7 2 ~ 2 1} j d cos 7 , (9) 

where we set the free parameter 7 m (the width 
of the distribution of directions) empirically. We 
repeat this 10 times for each history, resulting 
in a total of 5000 realizations. Unless specified 
otherwise we consider a halo at z = of mass 
5 x 10 n M Q in a ACDM cosmology with Q m = 0.3, 
n A = 0.7, H = 70 kms-^pc -1 , and a s = 1.0. 

3.2. Results 

The top-left panel of Figure 2 shows the result- 
ing distribution of A' for the choice / = 0.34 and a 
uniform distribution of directions 9. We note im- 
mediately that the distribution of A' is well fit by a 
log-normal distribution, equation 3, with values of 



Ag and <j\ fairly close to those seen in N-body sim- 
ulations, shown for reference in the bottom-right 
panel. 

We note that the value of o\ obtained with a 
fixed / and a uniform distribution of directions is 
somewhat larger then that seen in the simulations. 
This can be cured if we relax the assumption of a 
uniform distribution of directions, and allow for 
some correlation between the directions of succes- 
sive mergers, as detected in N-body simulations 
on scales of a few hundreds of kpc (Dekel et al. 
2000). We approximate this correlation by draw- 
ing the direction angle 9 from a Gaussian (eqn. 9, 
above). As seen in the upper- right panel of Fig. 
2, the choice j m = 2.0 reduces the value of <j\. 
The desired value of Ag and <j\ are now obtained 
by choosing a somewhat higher value of / = 0.38 
(Fig. 2, bottom-left panel). 

In reality, the value of / will also vary from 
encounter to encounter, resulting in a distribution 
with some scatter. This would tend to increase the 
scatter in A', which can be compensated for by a 
somewhat stronger correlation between the direc- 
tions, namely a further decrease in the value of 
7,„. For example if / is normally distributed with 
a standard deviation of 0.1, then to match the A' 
distribution requires j m = 1.0 and the mean value 
of / to be reduced to 0.32. Here, as we don't know 
the form of the distribution of / we will keep things 
simple by assuming it to have one value. Later, 
our recipe could be refined by extracting from N- 
body simulations the detailed distributions of / 
and 9, and their possible dependence on mass and 
redshift. 

4. The Tidal Torque Scenario 

The standard picture of the origin of the an- 
gular momentum of DM halos, as originally intro- 
duced by Peebles (1969) and then further devel- 
oped by Doroshkevich (1970) and White (1984), 
is formulated in the context of linear tidal torque 
theory. Here, material acquires its angular mo- 
mentum from the large-scale tidal field while still 
in the linear regime. We now develop an algorithm 
for the generation of angular momentum in virial- 
ized halos based on the tidal torque picture and a 
spherical collapse model. 

Linear theory (e.g. White 1984) predicts that 
the angular momentum of material before turn- 
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around, at time t, is given by: 

J i {t)=a(t) 2 D{t)e ijk T jl I lk , (10) 

where the time growth is from some fiducial initial 
time ij, Iik is the inertia tensor of the protohalo 
at fj, and Tji is the tidal (or shear) tensor at the 
halo center, smoothed on the halo scale. This is 
based on assuming the Zel'dovich approximation 
for the velocities inside the protohalo, and a 2 nd - 
order Taylor expansion of the potential. Porciani 
& Dekel (2001) have shown that the implied stan- 
dard scaling relation should be slightly modified. 
When applied to a collapsed shell of mass m and 
comoving radius q, it reads 

dJ = ga 2 (t c )D(t c )a(M h )mq 2 (11) 

where a(Mh) is the rms density fluctuation on 
scale of the mass interior to the shell at ij. The 
time t c is when the shell practically stopped gain- 
ing spin. An empirical fit from simulations shows 
that the effective time is, on average, about one 
third of the spherical- model collapse time, namely 
slightly before turn-around (Porciani, Dekel, & 
Hoffman 2001a). The parameter g represents the 
small mis-alignment between the inertia and tidal 
tensors and we expect it to be ~ 0.1 (Porciani 
et al. 2001b). 

We incorporate this into the merger trees by 
identifying the mass in a collapsing shell with the 
total mass of merging satellites (progenitor halos) 
and accreted mass during the same period of time. 
Of course the value of g will vary from halo to halo, 
so we assume that g is distributed normally with 
a mean go and a standard deviation a g . Based on 
the smoothness of the tidal field across the proto- 
halo volume, we assume that the angular momen- 
tum of each shell is in the same direction and can 
be added linearly. Since N-body simulations find 
that there is some nonlinear evolution of the spin 
direction at late epochs (Porciani et al. 2001b), we 
anticipate that g will be forced to vary with z. 

We again produce 5000 values of A' by consider- 
ing 500 mass accretion histories with 10 randomly 
chosen values of g. At z — in a ACDM cosmol- 
ogy, and for a halo mass of 5 x 10 u Mq, we find 
that the values go = 0.125 and <j g = 0AAg o yield a 
good fit to the A' distribution in simulations (log 
normal with X' Q = 0.35 and <j\ = 0.50). 
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Fig. 3. — The mean and scatter of the spin at 
z = as a function of halo mass, as derived from 
the two algorithms. The upper panel shows Aq 
(large connected symbols) and the corresponding 
Ao (small symbols) for the orbital-merger picture 
with / = 0.38 and j m — 2.0. The middle panel 
is the analog for the tidal-torque picture with 
go = 0.125 and a g — 0.44<?o- The bottom panel 
shows a\ for the two scenarios. The horizontal 
dot-dashed lines show the mass-independent val- 
ues found in Bullock et al. (2001), with the corre- 
sponding errors as error bars. The dotted line is 
the mean value of A from Bullock et al. . 

5. Independence of halo mass, cosmology 
and redshift 

An interesting property of the distribution of 
the halo spin parameter, A, as measured in cos- 
mological N-body simulations, is that it seems to 
be independent of halo mass, redshift and cosmol- 
ogy (Lemson & Kauffmann 1999). Any model of 
angular momentum acquisition should therefore 
produce the same spin distribution irrespective of 
these three parameters. Indeed, our simple mod- 
els seems to show this ubiquity of the A distribu- 
tion of halos, or can be adjusted to reproduce this 
independence by simple, physically motivated ad- 
justments of the free parameters. 

5.1. Halo Mass 

Figure 3 shows the values of X' Q as a function of 
halo mass with the orbital-merger picture (upper 
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panel, large diamonds) and the tidal-torque sce- 
nario (middle panel, large squares). Also a\ as a 
function of halo mass is shown for both scenarios 
(bottom panel). We see that over the mass range 
most relevant for spiral galaxies (10 10 — 10 12 M Q ) 
the distributions are all consistent with each other 
within the estimated errors (and the differences 
between different simulations). 

For the orbital-merger model there is a slight 
apparent trend for X' to increase and for o\ to de- 
crease as the mass increases. Note that when mov- 
ing from A' to A this trend is expected to weaken 
because the halo concentration parameter is also 
a function of halo mass. Using the trend of av- 
erage c v i r with halo mass found in (Bullock et al. 
2000a) we can convert each A' value to a value of 
A, though this ignores the spread found in values 
of the c v ; r and any possible correlations between 
A' and c v ; r . The corresponding average values of A 
are also plotted (upper panel, small diamonds) in 
Fig. 3 and show the absence of a trend with mass. 
In the simulations, there are some reports of an 
apparent trend of A with halo mass, in the oppo- 
site direction (Barnes & Efstathiou 1987), but this 
trend is weak and not seen in more recent simula- 
tions. 

For the tidal-torque method different mass ha- 
los as seen in Figure 3 are all in rough agreement 
with the N-body results, though there is a trend 
towards a decrease in X' and <j\ as mass increases. 
The trend in \' Q is in the opposite direction to 
that found in the orbital-merger scenario discussed 
above (which means that the dependence of con- 
centration on mass implies that the trend for A is 
even stronger (middle panel, small squares)). 

Thus both scenarios produce values of A' con- 
sistent with the results of simulations, though 
both seem to show some weak mass dependence. 
Whether this weak mass dependence can also be 
found in simulations remains to be seen. Since the 
trends with mass of the two scenarios are in op- 
posite directions, detection of such a trend would 
tend to favor one scenario over the other. 

Both trends can be removed by scaling one of 
the free parameters with mass. The functional 
forms 

f(M h ,z)= 0.35 (^J (12) 



go(M h ,z)= 0.17 , (13) 

where M* is the characteristic nonlinear mass (see 
Lacey & Cole 1993), create A' distributions with 
no mass dependence (note that a g remains 0.44g )- 
Our reason for scaling these expressions in terms of 
M* , which is a function of redshift and cosmology, 
will become apparent below (especially §5.3). 

For the orbital merger picture the weak correc- 
tion would imply that the incoming angular mo- 
mentum of satellites decreases slightly as a the 
halo's mass increases. In the tidal-torque picture 
we find that to remove the trend of A' with mass 
requires that the misalignment between the tidal 
tensor and the protohalo inertia tensor is larger 
for more massive halos. Both of these conditions 
can be checked in N-body simulations. 

We also observe a slight trend of o\ to decrease 
with increasing mass in both scenarios. This is a 
natural outcome of our modeling where the spread 
in A' values arises from the spread in halo forma- 
tion histories. It is unlikely for a halo to have a 
progenitor at a given redshift with a mass much 
greater then the value of M* at that redshift. Thus 
the more massive a halo is the less likely it is to 
have assembled most of its mass long ago. This 
naturally leads to a narrower range of possible his- 
tories for more massive halos, and thus to the ob- 
served trend of <j\ with halo mass. The effect is 
weaker in the orbital-merger algorithm where the 
spread in A' is also caused by the random orienta- 
tions of the incoming satellites. Future studies of 
simulations will be able to reveal whether such a 
trend exists. 

5.2. Cosmology 

To test the dependence on cosmology, we gen- 
erate merger trees based on two other cosmo- 
logical models: the original SCDM model, with 
an Einstcin-dcSittcr cosmology (f2 m = 1.0,51a = 
0.0, H = 50 kms^Mpc -1 ,^ = 0.5), and an 
OCDM model, with an open cosmology (Q m = 
0.3, ft A = 0.0, Ho = 70 kms _1 Mpc _1 ,CT 8 = 1.0). 
We compute the predicted A' distribution in each 
of these cosmological models using the same fidu- 
cial values of the parameters /, 7 TO , go and a g 
that produced the best fit to the A' distribution 
in the ACDM cosmology. The predictions of the 
orbital-merger scenario are found to be consistent 
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Fig. 4.— The distribution of A' in SCDM 
and OCDM cosmologies predicted by the orbital- 
merger (OM) and tidal-torque (TT) models. 

with the distribution of A' in simulations both for 
SCDM and OCDM. The tidal-torque scenario re- 
covers the distribution for the SCDM cosmology, 
but it predicts a somewhat higher value of X' for 
OCDM. These distributions slightly improve if the 
scaling with M* proposed in the previous section 
is used. Of course the distribution of A' in other 
cosmologies can be made identical to that of the 
ACDM cosmology by choosing slightly different 
free parameters. Again it will require higher res- 
olution simulations to see if such trends exist or 
not. 

5.3. Redshift 

Figure 5 shows the A' distribution as predicted 
by the two scenarios across the redshift range to 
3. We use the ACDM cosmology and the fiducial 
fixed values for the parameters /, j m , go and a g . 

In the orbital-merger case there is a slight in- 
crease of the mean value of A' with redshift. As 
in the case of differing masses, the change in the 
average value of the standard spin parameter A 
is weaker because concentrations evolve with red- 
shift as (1 + z)^ 1 . The mean value of A is also 
shown in Fig. 5, and shows almost no trend with 
redshift. It will be interesting to see whether the 
trend in the average value of A' is verified in N- 
body simulations. 
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Fig. 5. — The mean and standard deviation of 
the spin parameters for halos of mass 5 x lO n M0 
as a function of redshift, as derived by the two 
algorithms. The notation is as in Figure 3. Note 
that by redshift 1.5 concentrations become so low 
that A values overlap the A' values. 

The tidal-torque method with fixed parameters 
produces a noticeable evolution in the mean value 
of A', dominated by the fact that shells that col- 
lapse later have had more time to gain spin via 
tidal torques. The evolution is even stronger when 
A is considered because of the evolution of c V h- with 
redshift mentioned above. In order to keep the 
predicted A' distribution constant with redshift, 
one needs to vary the parameter go- Using the 
scaling of the free parameters / and go introduced 
above in terms of M^/M* (equations 12 and 13), 
the redshift dependence of the A' distribution is 
removed. 

In both scenarios we see a trend for o\ to de- 
crease with increasing redshift. Considering the 
mass of the halo over M* it is apparent that this is 
the same effect as seen with variation in halo mass; 
halos with masses closer to or greater then M* 
have less of a spread in formation history which 
leads to a smaller spread in their A' distribution. 

6. Merger History Dependence 

Having shown that the A distribution of halos 
resulting from our schemes can be made indepen- 
dent of cosmology, redshift or mass, in agreement 
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with what is found in N-body simulations, we now 
turn to where our scheme does predict a A depen- 
dence — the merger history. Preliminary results 
from N-body simulations suggest (Gardner 2000; 
Wechsler et al. 2001a) that the distribution of A is 
different for halos that have recently undergone a 
major merger. Using the same definition for a ma- 
jor merger (progenitor mass ratio of 1:3), we see 
a similar effect in our models. Figure 6 shows a 
significant difference between the distribution and 
mean value of A' for halos that have undergone a 
major merger since z = 0.5 in both scenarios. In 
work in progress based on a ACDM simulation, 
Wechsler et al. find a mean value of log A' = 0.044 
for halos identified at z = that have had a ma- 
jor merger at z < 0.5. Those that have not had a 
major merger since z = 0.5 have a log average of 
0.032, and those that have not had such a merger 
since z = 2 have a log average of 0.030. Gardner 
finds a similar result, a log average lambda value of 
0.049 for halos that have experienced major merg- 
ers since z = 0.5 in an OCDM model. 2 

Our two different scenarios predict similar val- 
ues for the mean in the log of the A' distribution for 
halos with recent major mergers, 0.062 and 0.048 
for the orbital-merger and tidal-torque scenarios 
respectively. The tidal-torque value is in good 
agreement with the simulations while the orbital- 
merger scenario gives a result « 30% higher then 
that found in simulations. If there is a situation 
where this algorithm might be suspect it would 
be the case of nearly equal mass mergers with 
high orbital angular momentum. In this situa- 
tion our assumption that we can neglect mass loss 
is probably invalid, and if mass loss occurs there 
will inevitably also be associated loss of angular 
momentum. However, recent major mergers are 
commonly held to result in ellipticals so the effect 
of our overestimating A' for these halos will have 
a negligible effect on disk galaxies. Our approach 
could be refined by treating the expected mass and 
angular momentum loss in a scheme similar to the 



A direct comparison between this simulation and our 
method is difficult because of differences in time step, halo 
identification, the concentration of halos, and the cosmo- 
logical model. These factors make it unsurprising that we 
find 1120 out of 5000 halos to be recent major mergers while 
Gardner finds only 15 out of 1609. Our results for the frac- 
tion of halos with recent mergers are in general agreement 
with the simulations analyzed by Wechsler et al. (2001a), 
where a more consistent comparison is possible. 



OM z<0.5 




OM z>0.5 



0.15 



Fig. 6. — Correlation of spin and merger history. 
The distribution of A' for halos that have and have 
not undergone major mergers since z=0.5 in both 
the orbital-merger (OM, top panels) and tidal- 
torque (TT, bottom panels) pictures is shown. It 
is surprising how similar the distributions are in 
these rather different pictures. We note that the 
log-normal distribution is not as good a fit to those 
halos that are recent major-merger remnants. 

one we have proposed. 

Mergers seem to leave a discernible mark on 
the halo's spin all the way back to z — 2. Ta- 
ble 1 shows the values of X' and <j\ for the two 
scenarios for halos that have and have not under- 
gone major mergers since redshifts of 0.5, 1.0 and 
2.0. 3 This strong connection between the merger 
and mass accretion history of a dark matter halo 
and its spin clearly has many interesting implica- 
tions for galaxy formation models. This will be 
best investigated in detail using full semi-analytic 
models, but we discuss qualitatively some of the 
implications below. 

7. Discussion 

We have presented algorithms for tracing the 
acquisition of spin angular momentum by dark 
matter halos through their mass accretion or 



3 Like Wechsler et al. and unlike Gardner, we find that the A' 
of non-mergers is different than the population as a whole. 
This is likely because we identify a significant fraction of 
halos as having undergone major merger events. 
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Mergers Nonmergers 
orbital merger tidal torque orbital merger tidal torque 





fraction of halos 


K 


o\ 


K 


ox 


fraction of halos 


K 


<J\ 


K 


o\ 


z = 0.0 












1.0 


0.035 


0.50 


0.035 


0.50 


z < 0.5 


0.22 


0.062 


0.27 


0.048 


0.47 


0.78 


0.031 


0.46 


0.032 


0.45 


z < 1.0 


0.43 


0.053 


0.33 


0.044 


0.47 


0.57 


0.027 


0.41 


0.032 


0.45 


z < 2.0 


0.68 


0.044 


0.41 


0.039 


0.48 


0.32 


0.023 


0.41 


0.028 


0.44 



Note. — The table lists value of \' and u\ for those halos that have suffered a major merger (3:1) since redshifts 0.5, 1.0 and 
2.0 and those that have not. Also listed are the number of halos that fall into each category. 



merger history. We have demonstrated that both 
methods reproduce the distribution of spin param- 
eter found in N-body simulations when a small 
number of physically motivated free parameters 
are tuned appropriately. 

We have proposed two such schemes, one based 
on the transfer of the orbital angular momentum 
of merging satellites to the internal spin of the 
halo, and another based on tidal-torque theory ap- 
plied to spherical shells of collapsing material. It 
is interesting that these two limits of the linear 
and non-linear regimes produce similar results. In 
principle, the tidal field around a halo should in- 
fluence both the halo's accretion history and the 
orbital angular momentum of satellites, so that in 
the orbital-merger model, the tidal field is still ul- 
timately responsible for the generation of angular 
momentum. If we knew the details of the den- 
sity field around the halo then we could calculate 
the orientation and magnitude of the angular mo- 
menta of the merging satellites instead of draw- 
ing them randomly in a somewhat arbitrary fash- 
ion, which is of course what an N-body simulation 
does. 

We obtain spin distributions with respect to 
halo mass, redshift and differing cosmologies that 
are consistent within the uncertainties found in 
N-body simulations. We also find weak trends of 
the A' distribution with mass and redshift that 
are in the opposite directions for the two scenar- 
ios. The dependence of concentration on mass and 
redshift found by Bullock et al. (2000a) creates a 
A distribution without a mass or redshift trend in 
the orbital- merger scenario, but increases the sys- 
tematic variation of A' with redshift for the tidal- 
torque picture. The discovery of such a trend in 



simulations would be a strong conformation of our 
modeling and its direction could be used to dis- 
criminate between the two scenarios. 

We also find a trend for o\ to decrease as mass 
or redshift increases. This can be understood in 
terms of the spread in formation histories of the 
dark matter halo. The more massive a halo is in 
terms of M^/M* the less likely it is to have been 
in place for a long time as that would require pro- 
genitors with masses larger then M* of the cor- 
responding time. The smaller spread in forma- 
tion histories in our modeling translates directly 
into a smaller spread in values of A'. This also 
explains the trend of A with mass and redshift 
found in the orbital-merger picture. As M/,/M* 
increases the prevalence of recent major mergers 
increases which as we have shown in §6 correlates 
with higher spin values. 

This effect should also be seen in the tidal- 
torque method; however, there are two other ef- 
fects in this case that go in the opposite direction. 
Recall from equation 10 that the angular momen- 
tum gained by each shell is proportional to a{Mh). 
For larger masses cr(M^) is smaller and this leads 
to lower spin values. Likewise, at higher redshifts 
shells have had less time to gain angular momen- 
tum. These effects dominate over the expected 
increase in A' because of late time accretion. 

Because we are unsure if such trends exist in 
N-body simulations, we also include simple scaling 
relations in M^/M* for one of the free parameters 
in each scenario (/ or go) that effectively removes 
the trends with mass and redshift. N-body simu- 
lations can be used to test whether g and / are 
functions of M/j/M*. 

We have obtained the interesting prediction, 
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supported by preliminary recent work with N- 
body simulations, that there is a strong correla- 
tion between spin value and merger history. In 
particular, the distribution of spins for halos that 
have experienced recent major mergers is skewed 
towards higher values. This effect is discernible 
for major mergers at least as far back in time as 
z <~ 2. Further study of N-body simulations will 
test how accurately our algorithms trace this cor- 
relation between spin and merger history (Wech- 
sler et al. 2001a). 

Vitvitska et al. (2001) use N-body simulations 
to measure the distribution of the free parameter 
we term / and find its most common value to be 
i=a 0.5 — 0.6. While this is higher then the value 
of / we have used (0.38), they also find that to fit 
the A distribution they need to assume that 25% of 
the angular momentum is lost during the merger 
event, implying that the effective value of / is very 
close to what we have found. Closer study of sim- 
ulations will be required before we understand the 
details of how orbital angular momentum is trans- 
fered to spin in mergers and what dependencies 
are involved; however, the simple prescription we 
have provided here should be adequate for many 
applications. 

These recipes can now be incorporated in semi- 
analytic methods that attempt to follow all the 
relevant physical processes of gas cooling and col- 
lapse, star formation, supernovae feedback, etc, 
over the formation history of galaxies. The self- 
consistent treatment of angular momentum acqui- 
sition that we have proposed is a significant im- 
provement over previous work, in which spins were 
assigned to halos at random, from the overall dis- 
tribution of A irrespective of the halo's merger his- 
tory. 

We can already anticipate some of the conse- 
quences of the predicted strong dependence of spin 
on merger history. As an immediate example, this 
may explain the small spread in disk sizes found 
by de Jong & Lacey (2000). They argue that the 
standard Fall & Efstathiou (1980) picture of disk 
formation leads to a relation Rd oc XLj between 
the disk size, the spin parameter and the near-IR 
luminosity, implying that the spread in disk size 
at a given luminosity is log-normal with standard 
deviation ~ a\. However, they deduce from an ob- 
served sample of late- type spirals o\ = 0.36±0.03, 
which is many "sigma" away from the N-body re- 



sult of ~ 0.5. However, Table 1 demonstrates 
that considering subsamples of the halo popula- 
tion, divided by merger history, results in lower 
values of a\. If we assume that late-type galaxies 
inhabit halos that have not experienced a major 
merger since a redshift of 1 or 2, then we find 
0.41 < a\ < 0.45, closer to the observed value. 
Thus, perhaps the smaller-than-prcdicted spread 
in observed disk sizes may be partially understood 
as the result of the special (quiescent) merger his- 
tory of dark matter halos that harbor late-type 
galaxies. The distribution could be further nar- 
rowed by assuming that halos with very small 
spins also form early-type galaxies or spheroids 
via instability processes (Mo et al. 1998; van den 
Bosch 1998). 

A second consequence of connecting Hubble 
type with merger history in this way is that the 
mean value of A' for the galaxies with a quiescent 
history (no recent major merger) is predicted to be 
10 — 15% lower then the mean of all halos. This 
implies a similar reduction in the prediction for 
the mean value of the disk scale length at a given 
luminosity. Whether such a reduction is problem- 
atic for the model will depend on the detailed ar- 
guments used to connect disk size to A, and will 
be explored in a future work. 
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